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ABSTRACT 


We demonstrate that a nonthermal distribution of particles described by 
a kappa distribution can be accurately approximated by a weighted sum of 
Maxwell-Boltzmann distributions. We apply this method to modeling collision 
processes in kappa-distribution plasmas, with a particular focus on atomic pro¬ 
cesses important for solar physics. The relevant collision process rate coefficients 
are generated by summing appropriately weighted Maxwellian rate coefficients. 
This method reproduces the rate coefficients for a kappa distribution to an esti¬ 
mated accuracy of better than 3%. This is equal to or better than the accuracy of 
rate coefficients generated using “reverse engineering” methods, which attempt to 
extract the needed cross sections from the published Maxwellian rate coefficient 
data and then reconvolve the extracted cross sections with the desired kappa 
distribution. Our approach of summing Maxwellian rate coefficients is easy to 
implement using existing spectral analysis software. Moreover, the weights in the 
sum of the Maxwell-Boltzmann distribution rate coefficients can be found for any 
value of the parameter k, thereby enabling one to model plasmas with a time- 
varying K. Tabulated Maxwellian htting parameters are given for specihc values 
of K from 1.7 to 100. We also provide polynomial hts to these parameters over 
this entire range. Several applications of our technique are presented, including 
the plasma equilibrium charge state distribution (CSD), predicting line ratios, 
modeling the influence of electron impact multiple ionization on the equilibrium 
CSD of kappa-distribution plasmas, and calculating the time-varying CSD of 
plasmas during a solar flare. 


Subject headings: atomic data, atomic processes. Sun: flares, techniques: spectroscopic 
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Introduction 


Evidence for no n-Maxwellian par ticle distributions has been found in the E arth’s 


magnetosphe re (e.g., 


corona (e.g., 


2009 


Vasyliunas 


Cranmer 


Oka et ah 


19681) ■ the solar wind (e.g., 


20141) . and solar flares (e.g, 


Feldman et ah 


Seely et al 


1987 


1975), the solar 


Kasparova fc Karlickv 


20131) . These nonthermal distributions are commonly characterized in terms 


of kappa distributions, which resemble a Maxwell-Boltzmann distribution at low energy but 
fall off as a power law at high energies. Kappa distributions are predicted by some statistical 
mechanical theories to be a natural consequence of systems in which the re is on-going 


heating, such as due to r econnection, shocks, or wave-particle interactions fjPierrard &: Lazer 


2010 


Dudfk et al 


2015 


and references therein). Thus, kappa distributions are likely to 


exist in a wide variety of astrophysical systems, even beyond the solar and space plasmas 
that have been the focus of existing work. 

The detection of nonthermal electron energy distributions (EEDs) via spectroscopy 
would be a powerful diagnostic. In solar physics, it has been argued that the presence of 


high energy electrons supports nanoflare theories of coronal heating flTesta et al 


29141 ) . 


Measurements that characterize the FED could indicate where such nanoflare heating 


is occuring. Such measurements may also pro vide insight in 


reconnection that lead to particle acceleration flRenz. fc Gudel 


o the processes wit hin 


29191) . 


Dudik et ah 


(129151) 


provide a review of recent results concerning kappa distributions in the solar corona. 

In order to detect the presence of kappa EEDs, or other non-thermal distributions, 
appropriate atomic data are needed to model spectra and analyze observations. For 
example, for an optically thin spectral line emitted by a tr ansition from level s to level / of 


charge state q for element X, the line intensity is given by (jPhillips et al. 


29981) 




( 1 ) 


















































4 


where is the electron density and the line element dh lies along the line of sight. The 
contribution function G^/ describes all of the atomic parameters for the transition and is 
dehned as 

n.iX"*) n(X'>+) n(X) «(H) A., 

n(X^+) n(X) n(H) Ue Ue 

Here, is the relative population of the upper level s for charge state X'^’*'. 

In collisionally excited plasmas, this level population is determined by the balance among 
collisional excitation, collisional de-excitation, radiative decays, cascades from the decays 
of higher energy levels, and recombination and ionization into and out of the various 
energy levels. The next term, n{X.'^~^)/n{X.), is the relative abundance for charge state q of 
element X. This describes the charge state distribution (CSD) of the plasma, which for a 
collisionally ionized plasma is determined by the balance among electron impact ionization 
(Eli) and electron-ion recombination. The other terms give the elemental abundance of X 
relative to hydrogen, n(X)/n(H), and the radiative transition transition rate Agf. 

A problem for interpreting spectra from plasmas with a kappa distribution of 
electrons is that the necessary atomic data have usually only been reported for Maxwellian 
distributions. The level population and the CSD both depend on the EED, as well as on the 
electron density. For a Maxwellian plasma, the influence of the distribution is characterized 
by the temperature T. However, for a kappa-distribution plasma, the level population and 
CSD also depend on the parameter k, which characterizes the degree of non-thermality of 
the distribution. 


There are several possible approaches available for obtaining the data needed to model 
collision processes in kappa-distribution plasmas. The most commonly used method has 
been to “reverse-engineer” the data that have been tabulated for Maxwell-Boltzmann 
EEDs so as to extract the required cross se ctions and then reconvolve these extracted cr oss 
sections with the desired kappa EED (e.g.. 


Dzifcakova 


1992 


Wannawichian et ah 


2003 
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Dzifcakova et al. 


20151) . An alternative method is to approximate the kappa distribntion 


itself as a snm of Maxwellians. With this latter approach, the needed atomic collision rate 
coefficients can be represented as a simple weig hted snm of the Maxwellian rate coefficients. 
This approach has been nsed, for example, by 


Ko et al. 


(1199611 to model the CSD in the 


solar wind. However, snmming Maxwellians has been relatively neglected in recent work. 


Here, we demonstrate that, in a general way, kappa distribntions can be approximated 
to a very high accnracy by a weighted snm of Maxwellians. We apply this general 
“Maxwellian decomposition” method to modeling collision processes in kappa-distribution 
plasmas. With this method, atomic process rate coefficients, generated by summing 
the appropriately weighted Maxwell-Boltzmann rate coefficients, can reproduce the rate 
coefficients for a kappa distribution to a level equal to or better than the accuracy obtained 
by using a reverse-engineering method. Moreover, the Maxwellian decomposition method is 
highly adaptable and can be readily implemented using standard plasma modeling codes. 

The rest of this paper is organized as follows: In Section |2] we introduce the kappa 
distribution. Section E] reviews the reverse-engineering approach to generating atomic data 
and then presents our methods and fitting parameters for representing kappa distributions 
as a sum of Maxwellians. In Section H] we compare our results for rate coefficients, CSDs, 
and predicted line intensities to those obtained using the reverse engineering approach. We 
then extend our analysis to several new applications, including ascertaining the influence of 
electron impact multiple ionization on the CSD for a plasma with a kappa distribution, in 
Section El and describing the time-dependent evolution of a plasma following a change in k, 
in Section [6l Section [7] concludes. 
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2. Kappa Distributions 


An energy distribntion quantifies the fraction of particles having an energy between E 
and E + dE. The isotropic Maxwellian distribution is given by 


f{E) = 




( 3 ) 


where k-Q is the Boltzmann constant and T is the temperature. 


The isotropic kappa energy distribution is given by 


f^{E) = 


TT V 


3/2 


\/E 


1 + 


E 


{k - 2>/2)kBT^_ 


— (k+1) 


with 


r(K + i) 

" r(/s:-l/2)(K-3/2)3/2- 


( 4 ) 


( 5 ) 


Here T is the Gamma function. The parameter k ranges from 3/2 to oo with k = oo 
being a Maxwellian distribution and smaller values of k corresponding to an increasingly 
non-thermal distribution. The kappa temperature is dehned so that the average energy 
of the particles is T^avg = (3/2)fcB^K5 in analogy with the usual Maxwellian temperature. 
The core of the kappa distribution can be approximated by a Maxwellian of temperature 
T = Tf^{K — 3/2) /k. In order for the magnitude of this Maxwellian function to match that 
of the ka ppa function at l ow energies, the Maxwellia n must be scaled by a multiplicative 


factor of flOka et ah 


2013 : 


Dzifcakova fc Dudfk 


20131 1 


ru + i) o/. 

C = 2.718—4- 

T K - 1/2 


1 \ 

1 + - 


( 6 ) 


This scaling factor can be thought of as representing the ratio of the number of thermal 
particles relative to the total number of particles. That is, it quantihes the fraction 
of particles at low energies where the distribution function more closely resembles a 
Maxwellian. Alternatively, Rn = 1 — G is the fraction of non-thermal particles, which is 
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largest for small k and decreases to zero for large k fjOka et al 


20131) . For a small value of 


K = 1.7 there is a large fraction of nonthermal particles with = 0.4, while a moderate 
K = 8 corresponds to = 0.1, and for a large value such as n = 100 the fraction of 
nonthermal particles is only R^ = 0.01. 


3. Atomic Data for Kappa Distributions 


In order to model the spectra from an electron-ionized plasma with a kappa FED it is 
necessary to determine the rate coefficients for the various atomic processes discussed above, 
in Section [1] These rate coefficients give the number of reactions per unit volume per unit 
time. As an illustration, the fractional ion abundance of charge state q is yq = n{X'^^)/nQCj. 
As a function of time, this is described by 

= rig Oj a^ q yq OR^g+l^/g+l] , (7) 

where of ^ is the rate coefficient for p-times ionization from charge state q to q + p and UR^g 
is the recombination rate coefficient from q to q — 1. The terms on the right of the equals 
sign represent, from left to right, ionization from lower charge states into q, ionization and 
recombination out of q to other charge states, and recombination from g -|- 1 into q. For 
most astrophysical atomic plasmas, the density is low so that three-body recombination 
is extremely unlikely. We also ignore charge exchange, which is only important at low 
temperatures of 10^ K. 


All of the needed ionization and recombination rate coefficients depend on the EED. 
Here and throughout we assume the distribution function to be isotropic. For any process 
having a cross section as a function of speed a{v) or energy (t{E) the rate coefficient a is 


given by (e.g.. 


Muller 


19991) 


I2E 


a = I a{v)vf{v)dv = / a{E)^l —f{E) dE, 


( 8 ) 









where fi is the reduced mass for the collision system. We are considering collisions of 
electrons with atoms and atomic ions so that /i ~ mg, although our general results also 
apply to other types of particle collisions. In the case of a Maxwellian distribution, the rate 
coefficients are generally a function of T. For kappa distributions, the rate coefficients are 
also a function of k. 

The problem for modeling plasmas with kappa EEDs is that atomic data are usually 
reported as Maxwellian plasma rate coefficients. If the cross sections themselves were 
given, it would be straightforward to perform the integral of Equation ([H]) to obtain the 
kappa-distribution rate coefficient. This would, in principal, be the most accurate way to 
obtain the needed rate coefficients. For some processes, such as electron impact ionization, 
the cross sections themselves are often available. However, for other processes, especially 
resonant processes, such as dielectronic recombination, it would be extremely cumbersome 
to tabulate the cross section in a format that could be incorporated into plasma models. 
Integration over a Maxwellian distribution smoothes out the complex structure in the 
cross-section data and allows for a simple parameterization. For this reason, it is Maxwellian 
atomic data that are usually tabulated, while keeping databases of cross sections is, for 
now, considered impractical. 


3.1. Reverse-Engineering Approach 


One method for obtaining rate coefficients for a nonthermal distribution is to 


extract the appropriate rate coefficients from the Maxwe 


fOwocki & Scnddei 

1983 

: Dzifcakova 

1992 

Hansen & Shlavantseva 

2004 

). The idea is 


Porqnet et al 


lian- i ntegrated atomic dat a 


2001 


Wannawichian et al 


2003 


rate coefficients for Maxwellian distributions and extract from them an approximate cross 
section, which can then be integrated over a kappa distribution. 
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To illustrate this reverse-engineering approach, we will discuss the method for 
dielectronic recombination (DR). DR is a process in which a free electron ap proaches an 


ion, excites a bound electron within the ion, and is simultaneously captured flMuller 


200811 


The resulting doubly excited state may relax by autoionization or by emitting a photon. 
Recombination occurs when the excited state relaxes radiatively to below the ionization 
threshold of the recombined system. DR occurs at collision energies E = AE — where 

AE is the electronic core excitation energy of the recombined ion and E\^ ~ eV 

is the bound-state energy of the captured electron into a Rydberg level with principal 
quantum number n of the ion with initial charge q. Because both AE and are quantized, 
DR is a resonant process. For each electronic state there are an infinite series of resonances 
corresponding to all possible Rydberg levels for the captured electron. Detailed calculations 
and mea surements of DR have been performed that resolve this complex resonant structure 


(see e.g.. 


Schippers et af 


2010 


, and references therein). 


The Maxwellian plasma rate coefficient aDR(T) is usually reported using the fitting 
formula 


«dr(7") — 




•I exp 


V^bT 


-E, 


(9) 


where bi and Ei are fitting parameters and the total number of terms in the fit is small, 
typically about ten. This approximation corresponds to a DR cross section that is a series 


of delta-function resonances at energies Ei, i.e., ctdr oc V., ctdr iS (E — Ei). 


coefficient for a kappa distribution can then be approximated as flDzifcakova 


'he DR rate 


1992h 


ttDR(7)c, k ) — 


A. 


(fcBT,)3/2 Z. 




1 + 


E, 


[k - 3/2)k 




— (k+I) 


( 10 ) 


Dzifcakoval fjl992h shows that this approximation is the first term of a series expansion and 


so the error can be estimated from the higher order terms. Qualitatively, the uncertainties 
are due to simplifying the complex resonant structure of the DR cross section into a 
small number of delta-function resonances. The largest uncertainties are expected for low 
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temperatures, where Ei ~ k^T and for small values of k —)■ 3/2. 


Similar methods have been given for approximating other releva nt processes, such 


Wannawichian et al. 

2003; 

Dzifcakova et al. 

2015 


Dzifcakova 


1992 


and de-excitation (e.g.. 


Dudik et ah 


2oii). 


2015!) and collisi o n stre ngths for excitation 


Dzifcakova et ah 


(12015!) recently developed 


a software package called KAPPA that tabulates kappa distribution rate coefficients for 
K = 2,3,4,5,7,10,15,25, and 33. This package integrates with the CHIANTI atomic 


database, w 

rich i 

Landi et al. 

2013) 


rich is widely used in the analysis of ultraviolet spectra (jPere et ah 


19971 : 


20131) . Below, we will compare our results using a Maxwellian decomposition 


approach to results obtained using the KAPPA package. 

The reverse engineering approach has several limitations. First, an adequate 
approximation must be found for each kind of atomic data, including RR, DR, electron 
impact ionization, collisional excitation and de-excitation, and so on. Including new atomic 
processes or improved methods of tabulating Maxwellian rate coefficients requires the 
development of new approximations. Updates to the Maxwellian data itself can require 
lengthy recalculations of the kappa distribution rate coefficients. For the case of DR 
discussed above, the changes are relatively simple substitutions, but for ionization or 
collisional excitation and de-excitation rates numerical integrations must be performed. 
Additionally , the a pproximations can become inaccurate for small k. For example. 


Dudfk et al 


fl2014fl discuss the approximation for distribution average collision strengths, in 


which the error in the approximation can be 20-30% for k = 2 and grows for smaller values 
of K. Finally, In order to model spectra using the rate coefficients derived in this way, an 
extensive database of all the required atomic data for each value of k of interest must be 
constructed. 
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3.2. Maxwellian Decomposition Approach 


An alternative approach to deriving atomic data for non-thermal distributions, is to 
approximate the distribution function itself as a sum of Maxwellians. We refer to this 


method as th e Maxwe 
For example, iKo et ah 


l ian de composition method. This method has been used in the past. 

axwellians 


(1199611 approximated kappa distributions using a sum of A 


in order to model the CSD in the corona and the solar wind. 


Kaastra et ah 


(1200911 


approximated the distributions of shocks in galaxy clusters using a sum of Maxwellians. 


If a distribution function can be approximated by a sum of Maxwellians, then the rate 
coefficient for the non-thermal distribution is just a sum of the Maxwellian rate coefficients. 
Suppose that the arbitrary distribution function g{E) is given by 


g{E)=Y,Cjf(E-,T,), (11) 

3 

where f{E]Tj) is a Maxwellian at temperature Tj, and the cj are parameters that weight 
the sum of Maxwellians. From the approximation of Equation (ITTl) and the expression for 
the rate coefficient given by Equation ([8]), it follows that the rate coefficient for g{E) is 

“g = (12) 

3 

where C({Tj) is the Maxwellian rate coefficient at the temperature Tj. Note that there is no 
requirement that the Cj be positive. 

The Maxwellian decomposition approximation has several potential advantages over 
the reverse-engineering approach. It is only necessary to calculate the approximation once 
for a given distribution function. Once the Cj and Tj are known, they can be applied to 
every atomic rate coefficient in the same way. Thus, the system is simple, extendable to any 
atomic process that can be represented as a Maxwellian rate coefficient, and does not need 
to be updated when new atomic data become available. The approximation for the needed 
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atomic data is, essentially, as good as the approximation of the Maxwellian decomposition of 
the kappa distribntion. Thus, with a suitable decomposition, data for small values of k can 
be obtained with accuracy as good as that for large values of k. Finally, it is straightforward 
to integrate into existing spectroscopic modeling software. In fact, the CHIANTI databas e 


already includes the ability to model spectra for a sum of Maxwellians (iLandi et ahl 120061) 


3.3. Method for Finding Kappa Fitting Coefficients 


In order to find fitting parameters that approximate a kappa distribution as a sum 
of Maxwellians, we have used two numerical procedures. We approximate the kappa 
distribution using a fit of the form 

f^{E] K, T^) = ^ Cjf{E- ajT^), (13) 

3 

where f{E]ajTi^) is the Maxwellian energy distribution at a temperature of T = a^T^. The 
Qj are independent of T^, but do depend on k itself. In order to further constrain the Cj we 
impose the normalization condition that 


= (14) 

3 

This ensures that the total integral of the distribution function is unity. 


One way to obtain the parameters Cj and aj is to simultaneously perform a least squares 
fit to the kappa function us ing a standard method, such as the Levenberg-Marquardt 


algorithm (jPress et ah 


19921). To determine the parameters, a set of energies Ei is chosen 


at which we want to match f^iEi). For the solutions given below, we used a linear, evenly 
spaced set of Ei starting at zero energy and extending up to an energy high enough so that 
less than 0.01% of the particles have an energy above the maximum Ei. In order to impose 
the normalization condition, all of the Cj are constrained to be greater than zero except for 
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one of them, which we label cat, which is set to 1 — order to provide the desired 

normalization constraint. We then performed the least sqnares £t to all of the parameters. 
We fonnd that this method works well for moderate to large valnes of k > 4. However, for 
more nonthermal distribntions, the hts tend to have a relatively large negative valne for 
Cat. As discussed below, in Section l3^ in order to reduce the uncertainties in the derived 
rate coefficients, the negative Cj should be kept small. Thus, we also used a different htting 
method, which minimizes any negative Cj values. 

This second method for performing the £t uses an iterative approach, in which we find 
the Cj and the aj separately. For a given set of aj, it is straightforward to hnd the best 
fit parameters Cj. Thus, we begin with an initial guess for the aj, solve for the Cj, then 
iteratively improve the aj and Cj as described below. The least squares best approximation 
to /k is the one that minimizes 


f^{Ei) - ^Cjf{Ei;ajT^ 


= Ti. 


(15) 


Here r* are the residuals between the Maxwellian decomposition and the kappa distribution. 
Since the aj are assumed to be known, and there are usually more energies than terms in the 
sum of Maxwellians, this is an overdetermined system of linear equations that can be solved 
using linear least squares, i.e., the solution minimizes The normalization condition 

of Equation fflTl) can be incorporated as a linear constraint. The best £t values for the Cj 


can then be found using standard methods (e.g.. 


Meyer 


19751) . We have used IDL, in which 


the laJeast-square-equality routine calculates the solution. There are equivalent functions 
in the linear algebra packages associated with many other programming languages. 


It is more difficult to optimize the set of aj. For our intial guess, we start with a few 
of the aj below one and the rest evenly spaced between one and the maximum E m^^ /knT^. 
Next the least squares solution for the Cj is found. For the purpose of minimizing the 
uncertainties (see Section l3^ . it is desireable that the magnitudes of the weights \cj\ be 
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small. This can be accomplished by minimizing the \cj\. 

We have fonnd that Ylj \^j\ is strongly correlated with the residnal error between 
/k and the Maxwellian decomposition. Fignre [1] shows an example of this correlation for 
K, = 7. Varions possible sets of aj were selected nsing a random nnmber generator and the 
corresponding Cj were fonnd. We then calcnlated the relative error 


Rc^xr — 




Ik 


and plotted the maximnm relative error for E < T^max 

Emax max -[ R^rr ) E < T/j^ax } 


(16) 


(17) 


against the \cj\. The Spearman rank correlation coefficient for this plot is 0.86 with at 
a significance of greater than 99.99%. Althongh we are not certain what the nnderlying 
reason for this correlation is, we can regard it as an empirical resnlt. 


Using this correlation, we can employ a minimization procednre t o find the set of a 


that minimizes Yj \^j\- do this nsing a downhill simplex method fjPress et al.lll992h . 


The resnlt does depend on the initial gness for the aj and in some cases even after the 
minimization procednre the resnlts still exceed the desired maximnm relative error. In this 
case, we can nse a random nnmber generator to pertnrb the initial gness and then repeat 
the procednre nntil we find a solntion that meets onr reqnirements. 


One other factor to consider, nsing either fitting method, is the nnmber of terms in the 
snm. We find that a larger nnmber of terms is needed for smaller valnes of k. Onr approach 
to this issne was to start with some moderate nnmber of terms, and then add or remove 
terms from the snm nntil it appeared that we had reached the minimnm nnmber needed for 
the desired accnracy. 

Fignre [2] shows an example of a fit for k = 2 at = 10. Note that the nnits of 

temperatnre are arbitrary as long as has the same nnits as E. The dashed line 
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in the figure shows and the solid line shows our Maxwellian decomposition 

approximation. The lower panel of the figure illustrates i^err, which is below 2.5% for 
E < 330 k^T^. This energy is indicated by the vertical dotted line on the plot. Fewer than 
0.01% of the particles in the distribution have energies greater than this. 


We have used both of the above fitting methods to find the parameters for the 
Maxwellian components, and report here the results that were best in terms of having a 
small error and \cj\. The procedure for finding a set of fitting parameters involves some 
trial and error. However, once a solution with sufficient precision is found it can be used for 
all of the atomic data; the procedure only ne eds to be performed once. A different numerical 


approach was used by 


Kaastra et ah 


(1200911 to approximate an arbitrary distribution 


function by a sum of 32 Maxwellians. The accuracy of their fits is comparable or sometimes 
worse than ours. Despite these issues, it is remarkable how well kappa distributions can be 
approximated by a sum of Maxwellians. We speculate that with a deeper understanding 
of the mathematical properties of kappa distributions a more systematic approach to this 
decomposition could be found. 


3.4. Uncertainties 

The uncertainties in the rate coefficients using the Maxwellian decomposition arise 
from the goodness of the approximation of the kappa distribution and from the weighting of 
the Maxwellian terms. This is in contrast to the uncertainties with the reverse engineering 
method, in which the uncertainties are due to the approximations involved with the 
individual rate coefficients. Both methods, of course, suffer from the same uncertainties in 
the Maxwellian atomic data. That is, if the tabulated Maxwellian data is inaccurate, then 
the derived data for the kappa distribution will also be inaccurate, though not necessarily 
in exactly the same ways. Since our interest is in comparing the decomposition and 
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reverse-engineering methods, we will omit discussion of errors in the Maxwellian data. 

One source of uncertainty for the Maxwellian decomposition is the accuracy of the 
approximation to the kappa distribution. This uncertainty can be mitigated by choosing a 
desired tolerance when finding the fit parameters. Below, in Section 13.51 we describe some 
cases where we require the maximum error to be less than 3%. In many cases we were able 
to obtain even smaller errors. 

Although we characterize the error in the Maxwellian decomposition using the metric 
Rmax, the error in the generated kappa-distribution rate coefficients can be different from 
this value. The error in the rate coefficient is a weighted average over the relevant cross 
section of the error at all energies, whereas the maximum Re„ occurs at the highest energies. 
If the cross section is largest at low energies and decreasing towards high energies, as is 
usually the case, the error in the rate coefficient is expected to be smaller than i?max- In 
some cases the cross section is at its maximum at very high energies, such as for some 
multiple ionization processes. In those cases, i?max niay underestimate the error in the rate 
coefficient. 

An uncertainty also arises due to the truncation of the fit at some maximum energy. 
That is, the fits are accurate to the desired tolerance for an energy range from = 0 to 
some maximum energy T'max- In our typical hts, we chose this energy to be such that no 
more than 0.01% of the particles have an energy greater than T^max- Because the cross 
sections for ionization, recombination, and excitation all decrease at sufficiently large 
energies, this truncation is not expected to produce signihcant errors in the generated 
relevant rate coefficients. 

A more complex source of uncertainty comes from the weighting of the terms in the 
sum. The Maxwellian rate coefficients are functions of temperature, and may have an 
uncertainty that is also a function of temperature, so that a Maxwellian rate coefficient 
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a{Tj) has an uncertainty of Sa{Tj). Assuming they are uncorrelated, the uncertainty in the 
derived kappa rate coefficient isar 

( 18 ) 

This implies that if any of the magnitudes of the Cj are greater than one, the uncertainty 
in the kappa-distribution rate coefficient must be greater than that of the Maxwellian rate 
coefficient from which it is derived, even if the kappa distribution itself is approximated to 
a very high accuracy. This is the main reason that we minimize \cj\ in our method. 
For our hts, none of the Cj is greater than one and the sum does not exceed unity by more 
than a few percent and usually by less than one percent. The resulting uncertainty in the 
kappa-distribution rate coefficient is then even smaller, because Iffi l- 

Another possible source of uncertainty is related to the accuracy of the tabulated 
Maxwellian rate coefficients. The decomposition method uses a sum of Maxwellians at 
various temperatures above and below to approximate the kappa distribution. However, 
the atomic data may be tabulated using a function that is expected to be valid within a 
certain temperature range; for CHIANTI this is typically 10^-10® K. The terms in the sum 
that fall at lower or higher temperatures may be inaccurate if the Maxwellian data are 
not accurate at those energies. Because the kappa distribution has a high energy tail, the 
decomposition is likely to have terms that greatly exceed T^, but usually it does not have 
terms at temperatures so much smaller than the minimum tabulated temperature as to be 
problematic. The errors at high temperatures are also expected to be small. This is because 
the Cj at those temperatures are small and because the Maxwellian rate coefficients are 
usually tabulated using a form that has the correct high energy behavior so that the errors 
when the rate coefficients are extrapolated to high temperatures are not grossly inaccurate. 
One caveat, though, is that the tabulated atomic data usually ignore relativistic effects, 
which may be important for very high temperatures. 
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Finally, because the Cj can be negative it is possible in some situations to find a small 
negative rate coefficient, which is unphysicah As shown in Section 1X51 the negative values 
of Cj are usually associated with S> 1. Negative rate coefficients tend to be produced for 
very low in processes that have a threshold, such as ionization. In this case, the largest 
values of Cj are concentrated at low temperatures where a^ajT^) is small, but a negative 
Cj can occur at a high temperature where the value of a{ajTf^) is large, resulting in a net 
negative rate coefficient. In practice, we have found these negative rate coefficients to be 
insignificant. They only occur in cases where the process is unimportant. For example, with 
ionization, negative rate coefficients occur at low where the abundance of the relevant 
ion is very small. Thus, any negative rate coefficients that arise can be set to zero with 
negligible consequences for the CSD or the level populations. In the unlikely case that 
the negative rate coefficient is problematic, the reader can use the method described in 
Sectionto find an improved Maxwellian decomposition with even smaller \cj\. 


3.5. Coefficients for Certain Values of k 

We have found fitting coefficients that describe kappa distributions as a sum of 
Maxwellians to very high accuracy. Tables HHSl list the aj and Cj parameters for k, = 1.7, 

2, 3, 4, 5, 7, 10, 15, 20, 25, 30, 33, 50, and 100. The coefficients given in these tables can 
also be obtained using the IDL code that is included as supplementary online material with 
this paper. The kappa distribution is approximated by substituting these parameters into 
Equation ([13]). 

TableElgives some parameters describing the accuracy of the Maxwellian decompositions 
in Tables [TH51 The accuracy is calculated for E < Emax, which is the energy below which 
99.99% of the particles in the kappa distribution are found. The third column of Table [6] 
lists the maximum error i?max- Also given is the sum of the magnitudes of the \cj\. These 




are all very close to unity with all of the cj < 1. This implies that the intrinsic errors in the 
Maxwellian rate coefficients will not be magnihed when calculating the kappa distribution 
rate coefficients. 


3.6. Coefficients for 1.7 < k < 100 

We have also found that it is possible to find an approximation in terms of a sum 
of Maxwellians that varies reasonably smoothly as a function of k. This allows for the 
possibility of varying k continuously. In order to find such fits, we fixed the set of aj values, 
based on what we found by optimizing the results for particular values of n (see Section 13.51 
above). We then performed the linear least squares analysis to determine the best £t Cj 
corresponding to that set of aj. It turns out that the Cj then vary smoothly as a function of 
K, and these Cj{K) can be fit by a polynomial. The resulting fits are accurate within a range 
of K values, so once the desired accuracy is no longer met, a new set of aj can be used and 
the results for the various Cj patched together piecewise. The end result is a continuous 
approximation of kappa distributions as a Maxwellian decomposition for k = 1.7-100. The 
results cannot be expected to be as accurate as the fits found by optimizing the aj for each 
K. However, the approach enables one to study how properties of spectra vary continuously 
as a function of k, though at the cost of some precision. 

Table [7] gives the aj and fitting parameters for the Cj in terms of polynomials in k, i.e. 

Cjin) = '^dnK"'. (19) 

n 

The aj and Cj reproduced using these polynomial fits can also be obtained by using the 
IDL code that is included as supplementary online material with this paper. When these 
and Cj are used to approximate the kappa distribution, the results have i?max ^ 25% in the 
energy range where 99.9% of the particles are found. The i?max is typically better than 30% 
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for the extended energy range that includes 99.99% of the particles. Figure |3] illustrates the 
accuracy as a function of k. Given that the uncertainties in the atomic data themselves are 
typically about 20%, this is a reasonable accuracy for modeling spectra. Additionally, the 
magnitudes of the \cj\ sum to less than 1.3 in all cases, and for most k the sum is below 
1.02. 

An alternative method for modeling spectra for continuously varying k is to perform 
a linear interpolation of the rate coefficients between two k values at which the rate 
coefficients are known accurately. Either method can be used, and the difference between 
the methods gives an estimate of the uncertainty in the modeled data. We have found small 
differences of < 10% for a few test cases of using linear interpolation of the rate coefficients 
versus the Maxwellian sum method with the Cj from the polynomial fits in Table [71 


4. Comparison of Methods 

As will be shown below, the Maxwellian decomposition and reverse-engineering 
methods give similar results for generating kappa-distribution atomic data and predicting 
properties of the plasma, such as the the CSD and line intensity ratios. Summing 
Maxwellians has, additionally, several other qualitative advantages. A major one is that 
it can be implemented in a straightforward way using existing spectroscopic analysis 
software, such as CHIANTI. This is in contrast to reverse-engineering methods, where 
a new database of rate coefficients must be constructed for each value of k. Also, with 
the decomposition approach, the parameters for the Maxwellian sum can be optimized 
until an approximation with the desired properties has been obtained. By adding more 
components to the sum, the approximation can be made accurate to very high energies. 
As the accuracy of the kappa distribution rate coefficients is mainly determined by the 
precision with which the Maxwellian sum approximates the kappa distribution, it should be 
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similar for all of the atomic data and any improvements in the decomposition approach are 
antomatically propagated into the derived data. This also implies that the approximation 
for very small valnes of k can be made as accnrate as that for larger, more Maxwellian, 
kappa distribntions. This is in contrast to the reverse-engineering approach, where the 
errors increase for small k, becanse some of the approximations nsed to generate become 
increasingly inaccnrate. Below, we compare onr resnlts for the rate coefficients, CSD, and 
emissivities with those obtained ns ing the reverse engineer ing approach. Specihcally, we nse 


the KAPPA package developed by 


Dzifcakova et al. 


fl2ni5h . 


4.1. Rate Coefficients 


KAPPA is a database of reverse-engineered rate co efficients and rela ted atomic data 


derived from those kappa-distribntion rate coefficients. 


Dzifcakova et al. 


have, for a few 


cases, compared of their rate coefficient resnlts with a direct integration of the cross sections. 
For those test cases they fonnd that the agreement was better than 10% and in most cases 
the precision is better than 5%. In order to give examples from both both relatively smooth 
and highly resonant processes, we consider ionization and recombination rate coefficients. 


For the ionization rate coefficients, we generally find excellent agreement between the 
reverse-engineering and Maxwellian decomposition approaches. FignreH] shows an example 
of the ionization rate coefficient for Fe^"*" forming Fe^"*" with k = 2. We choose k = 2 
becanse it is the most nonthermal distribntion in the KAPPA package. In the fignre onr 
resnlt is shown by the solid cnrve, while the dashed curve indicates the values given by 
the KAPPA package. The lower panel shows the relative error (the KAPPA-package result 
— our result) / our result. In this case we find that the derived rate coefficients agree to 
better than 5% at every temperature. We do, however, find that the KAPPA package gives 
a value that is systematically smaller than ours, by a few percent. Such discrepancies are 
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common for the other kappa rate coefficients and also occurs for the Maxwellian [k = oo) 
rate coefficient given in that package. For this reason, we suspect that the discrepancy 
is due to the truncation of the integral over the cross section in the KAPPA package 
calculations. Additionally, for some of our hts for different k, we hnd discrepancies at low 
Tfi,. This appears to be due to the negative Cj in our approximation. While these can lead 
to apparently large relative errors, we hnd that they have negligible inhuence on the CSD. 

The recombination rate coefficients also usually show good agreement, although 
there are some discrepancies at low for low charge states. Figure E] shows the total 
recombination rate coefficient (DR + RR) for Fe^"*" forming Fe^+ with k = 2. The solid 
curve illustrates our result, which is compared to the result from the KAPPA package as 
shown by the dashed curve. The lower panel gives the relative difference between the two 
methods. For this ion we hnd that there is a discrepancy of about 30% at low temperatures. 
A similar discrepancy is present for other values of k, which suggests that the error is 
caused by the approximations used in the reverse-engineering approach. Similar errors are 
found for other low charged Fe ions, up to about Fe®^, but for more highly charged ions 
the agreement is excellent, within a few percent, over the entire temperature range. For 
the Fe^"*" rate coefficient shown in the hgure, one can also see a discrepancy at very high 
temperatures. The reason for this discrepancy is unclear. However, as the abundance of 
Fe^^ is essentially zero at such high energies, this uncertainty will not affect the CSD. 


4.2. Equilibrium Charge State Distribution 


The ionization and recombination rate coefficients affect the spectra through the CSD. 
In order to further compare the Maxwellian decomposition and reverse-engineering, we have 


calculated the equilibrium CSD using our rate co efficients, ignoring mu 


the moment, and compare to the results given by 


Dzifcakova &: Dudfkl fl2013f) and available 


tiple ionization for 
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in the KAPPA package. 


The equilibrium CSD is known as collisional ionization equilibrium (CIE). Because the 
CSD is not evolving in time, the left side of Equation ([7]) is zero. This also implies that the 
density is a constant factor and plays no role in the solution. For a given temperature, we 
have a system of algebraic equations. It is easy to see that Equation ([7]) can be written as 
a matrix 

Ay = 0, (20) 

where A is the matrix of the rate coefficients and y is the vector of abundances with 
elements yq. Equation ([7]) incl udes multiple ionizat i on ra te coefficients, which were not 


included in the calculations of 


Dzifcakova &: Dudfkl (120131) . so for this comparison we will 


consider only single ionization. The effect of electron impact multiple ionization on the 
CSD for a kappa distribution is explored in more detail below, in Section [5l 

In order to obtain a unique solution to Equation fl20|) . an additional equation is needed. 
For this, we require that abundances be normalized so that 




( 21 ) 


This condition is implemented 
Equation fl2T|) . see for example 


by replacing one o f the rows of Equation (ITU]) with 


Brvans et ah 


(1200611 . 


For temperatures of = 10^-10® K, we find good agreement between our CSD results 
and those in the kappa package, with discrepancies of < 10%. However, at lower and 
higher temperatures there can be larger differences. The top panel of Figure [6] shows the 
relative abundances of Fe charge states for k = 2, with our results illustrated by the solid 
curves and those in the KAPPA package indicated by the dashed curv es. The lower pan el 


shows the ratio of our (New) calculations to the (Old) calculations of 


Dzifcakova et al. 


(|2015h . The curves are plotted only for those temperatures where the abundance of the 
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corresponding charge state is greater than 1%. Over most of the temperature range, the 
agreement is very good, which is also the case for other values of k. However, at very 
low and very high temperatures, the discrepancies can be ~ 50%. At low temperatures, 
one cause of the discrepancy is the inaccuracy in the reverse-engineering generated DR 
rate coefficients. The inaccuracy of the Maxwellian decomposition in the far tail of of the 
kappa distribution can lead to inaccuracy in the ionization rate coefficients and so can also 
contribute to discrepancies at low temperatures, as discussed in Section 13.41 The reason for 
the discrepancies in the CSD at high temperatures is not clear. 


4.3. Level Populations and Line Emissivities 


In addition to modeling the CSD, the Maxwellian decomposition method can also 
be used to model the level populations and emissivities of spectral lines. This modeling 
can be done in a straightforward way using the CHIANTI emiss_calc function with 
the sum_mwLcoeff keyword set. This function calculates the emissivities, that is level 
populations multiplied by radiative decay rates, for a given temperature and density. 
The keyword, causes the emissivities to be calculated for a di stribution that is a sum of 


Maxwellians weighted by the coefficients Cj (iLandi et ah 


2006 ). The KAPPA package also 


has a function for calculat ing emissivities t 


rat is an extension of the CHIANTI function 


and is called emiss_calc-k (IDzifcakova et ah 


20151). 


Here, we compare our results for the intensity ratio of several transitions of O vi, which 
is an abundant ion that is commonly observed in solar physics. The transitions we consider 
are ls^2s ^S' 1/2 3p ^Pi/ 2 , 3/2 at 150 A, ls^2p ^Pi/ 2 , 3/2 ^ 3d ^ 113 / 2 , 5/2 at 173 A, and 
ls^2s ^51/2 2s^2p ‘^Pz /2 at 1032 A. Figure [7] plots the ratios of the 150 A/1032 A and 

173 A/1032 a lines as a function of k for a fixed temperature of = 1 MK and density 
He = 10® cm“^, which are typical of the solar corona. The solid curves in the figure were 
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obtained using our method, with the Maxwellian distribution parameters found for every 
value oi K = 1.7-50. The hlled circles indicate the ratios given by the KAPPA package the 
values of k available in the package. 


Diagnostics for k, based on line intensity ratios have been proposed (IDzifcakova fc Kulinova 


2011; 

Dudfk et ah 

2014, 

2015) 


20151) . One limitation of the reverse-engineering method has been 


that it is a complicated process to obtain all of the atomic data necessary to calculate 
these ratios, and so only data for certain discrete values of k are available. The Maxwellian 
decomposition approach offers the possibility of easily modeling line intensity ratios for any 
value of K. Thus, given a set of observed line ratios, it is possible to hnd a parameter space 
of K, and Ue that is consistent with the data, and constrain all three of these quantities. 


5. Influence of Electron Impact Multiple Ionization on the CSD 


Electron impact multiple ionization (EIMI) is a process in which a si ngle electron-ion 


collision causes two or more bound electrons to be removed from the ion. 


Hahn fc Savin 


(120151) tabulated cross section data for EIMI of Fe, which can be used in CSD calculations. 


They showed that for Maxwellian plasmas, EIMI has only about a 5% affect on the CIE 
abundances, but for dynamic plasmas EIMI can signihcantly change the evolution of 
the plasma compared to what is expected if only single ionization is considered in the 
calculations. 


In plasmas with a kappa EED, the affect of EIMI on the CIE abundances may be much 
larger than for a Maxwellian plasma. This is because the kappa distribution has a long high 
energy tail. For a Maxwellian electron distribution, there are few electrons at energies above 
the energy threshold for EIMI. Thus, for a Maxwellian, very high temperatures are needed 
for EIMI to occur, but at such high temperatures the abundance of the ion undergoing 
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EIMI is very small. As a result, the net effect of EIMI on the equilibrium CSD is small for 
a Maxwellian distribution. In contrast, the high energy tail of kappa distributions implies 
that a signihcant fraction of nonthermal electrons have energies above the EIMI threshold, 
even at moderate temperatures. 

We have calculated the CIE abundances for iron including EIMI by solving Equation ([7]) 
with the time derivative set to zero, as described above in Section 14.21 However, in this 
case we keep the multiple ionization terms. The results for a large value of k = 25 and a 
very small k = 1.7 are shown in Figures [Eland [HI respectively. The solid curve in the figures 
shows the CIE abundances as a function of with EIMI included in the calculation, while 
the dashed curve illustrates the result if EIMI were ignored. The lower panel shows the 
ratio of the results that include EIMI to those that use only single ionization. In both plots, 
the curves are only plotted for where the abundances are greater than 1%. 


The figures show that EIMI has a signihcant effect on the CSD even in equilibrium for 
kappa distributions. Even for relatively weak kappa distributions, such as k = 25 the CSD 
differs by up to 40% from what would be expected with only single ionization. In this case, 
the largest relative effect appears to be in the low charge states, which become ionized to 
higher charge states at lower temperatures when EIMI is included. The inhuence of EIMI 
grows as the distribution becomes more nonthermal. For k = 1.7, the largest discrepancies 
occur for between about 10® and 10^ K, where EIMI can change the CSD by a factor 
of 2 to 7. The iron charge states in this temperature range have their outermost electrons 
in the M shell. EIMI signihcantly modihes the CSD at these temperatures, because the 
kappa distribution contains enough high energy electrons to ionize from the L shell. The 
resulting L -shell hole leads to a utoionization of an M-shell electron, ending in a net double 


ionization (iHahn fc Savin 


20151). 
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6. Transient Evolntion of a Plasma with Varying k 
6.1. Timescales 


The availability of rate coefficients for any value of k allows one to model plasmas with 
a changing k. One useful way to quantify the response of plasma to a change in k or is to 
calculate the timescale for the CSD to reach equilibrium following a change in the EED. For 
Maxwellian EEDs, these timescales are calculated for changes in T, and are useful in the 


Masai 

1984; 

Hughes & Helfand 

1985; 

Smith & Hughes 

2010) 


The same methods for calculating these timescales for changes in T in Maxwellian 
plasmas c an be u sed to derive the timescales for changes in and k in kappa plasmas 
Following iMasail fjl984l) . Equation ([7]) can be written as 

dy 


dt 


= n,A(K, T^)y. 


( 22 ) 


If K, Tk and He are constant, such as following a jump from a different set of values for these 
parameters, then the solution is y{t) = ^oexp [netA{K,Tf^)]. The exponential of a matrix 
is defined by the Taylor expansion in powers of the matrix A. The matrix multiplications 
are more easily carried out by diagonalizing the matrix by finding its eigenvalues and 
eigenvectors. Doing so, one finds that the solution to Equation (12^ can be written as 


y{t) = S 


exp (AiUet) 0 

0 exp (XkUet) 0 

: 0 


s m, 


(23) 


where S is the matrix in which each fcth column is the eigenvector corresponding to the 
eigenvalue of the diagonal matrix. The density-weighted timescales for equilibration are 
then given by 1/Xk in units of net. Note that the eigenvectors represent linear combinations 
of the elements of the vector y, rather than individual elements (i.e., charge states). 
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These timescales could be calculated for any combination of changes in k and T^. 
Here, we show the timescales for a change in k, while keeping hxed. Figure fill] shows 
the maximum timescale 1/Xk for as a function of k for = 10^, 10®, and 10^ K. This 
timescale represents the time it would take for the slowest evolving eigenvector to come to 
equilibrium. The hgure shows that the equilibration timescale is shorter for smaller values 
of K, where the distribution is most nonthermal, while the timescale is longer for large k 
where the distribution is closer to a Maxwellian. An explanation for this is that kappa 
distributions have both more low energy particles and more high energy particles than a 
Maxwellian with T = T^. This property is becomes more pronounced at smaller values of 
K. Since the recombination rates are generally greater for lower energies and the ionization 
rates tend to be greater for higher energies, distributions with more high- and low-energy 
particles respond faster. 


6.2. CSD Evolution for Solar Flares 


Nonthermal EEDs are produced in solar flares fe.g.. lSeelv et al.lll987l) . Some studies 
suggest t hat the energy distributions could be well des cribed by kappa distributions with 




2-8 (iKasparova fc Karlick 


ickJ 


2009 


Oka et ah 


3 ). 


20131) . Using our results for the atomic 


data for a kappa distribution, we can study how the CSD of the plasma evolves in response 
to a flare. 

In order to describe the CSD as a function of time, we solve Equation [7] numerically. 
Initially, we begin with a CSD that corresponds to a Maxwellian distribution, k = oo, 
with Tk = 1 MK. This approximately represents the quiescent state of the corona. We 
assume that the flare turns on as a sudden jump that immediately changes the EED 
and consequently the rate coefficients, but that t he CSD evolves more s l owly. For the 


parameters of the flare EED, we use results from 


Kasparova fc Karlickvl ( 2009 ). They 


















29 


analyzed data from hard X-ray spectra and found that for a coronal source, the spectra 
could be described using a kappa distribution with k ^ 8 and 15 MK. We used these 
values in our calculation. 


Figure [TT] shows the time evolution of several selected charge states of iron, following 
the sudden change in k and T^,. Three versions were considered. The solid curves show 
the case in which EIMI is included, k = 8, and = 15 MK. The dashed curve illustrates 
the same scenario, but includes only single ionization. Finally, the dotted curve includes 
EIMI, but models only the jump in and constant k = oo (i.e., the Maxwellian case). 
Comparing Erst the kappa curves with and without EIMI, it is clear that EIMI increases the 
rate at which the CSD evolves through the lower charge states, but has less of an effect for 
the higher charge states. Next, comparing the CSD evolution with k = 8 and EIMI versus 
the Maxwellian, we find several differences. First, the CSD evolves towards a different 
equilibrium, as can be seen in the plot at large times where the abundances of Fe^®"*" and 
Fe^^"*" are asymptoting towards different equilibrium values for the two distributions. The 
two distributions evolve on similar timescales, although in some cases the k = 8 CSD 
lags the Maxwellian CSD. This seems to contradict the timescale analysis above, which 
implied that the k distribution reaches equilibrium faster than a Maxwellian. However, the 
timescale analysis considers the general CSD rather than any particular charge state and 
so is not directly comparable to these results. For coronal densities of 10®-10® cm ^ the 
difference in the evolving timescales are of order seconds or tens of seconds. Such differences 
are at about the current tem poral resolution of solar observations, s uch as those from the 


Solar Dynamics Observatory flLemen et ah 


2012 


Woods et ah 


20121 ). 


Although we have focused on the parameters relevant to a coronal X-ray source, 
we have found qualitatively similar results for footpoint sources. The values of k at the 
footpoints that have been inferred from observations are usually smaller, with values of 
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K ~ 2-3. The analysis at the footpoints is also more complex as a model of the spectra 
should include the structure of the chromosphere and transition region, opacity, and 
hydrodynamic effects in more detail. Nevertheless, our results indicate that the evolution 
of the plasma CSD in response to a flare can be signihcantly changed by the presence of a 
non-thermal distribution, both in the corona or at the footpoints of the flare loop. 


7. Summary 


We have shown that a kappa distribution can be approximated by a sum of Maxwellians 
to an accuracy of better than 3%. This approximation can easily be used to generate 
collision data needed to model the CSD and other spectroscopic properties of plasmas with 
kappa EEDs. The resulting data are at least as accurate as those obtained using other 
methods. In particular, we have foun d that our results are in good agreement with those 


obtained using the KAPPA package of 


Dzifcakova et al, 


(120151) . This not only demonstrates 


that our Maxwellian decomposition method is a good alternative to the reverse-engineering 
approach, but also provides an independent estimate of the uncertainty in that method. 
For many purposes, either method can be used to a similar level of precision. 


However, summing Maxwellians has some advantages compared to other approaches. 
These are, hrst, that it is easy to use and can readily be incorporated into existing 
spectroscopic analysis software. Second, Maxwellian data for essentially any process can be 
used to generate kappa-distribution data by using the same small set of htting parameters. 
This contrasts with reverse engineering methods where specihc approximations must be 
used for each atomic process and large databases of rate coefficients must be generated 
for each k, in order to model spectra. A related advantage is that any updates in the 
tabulated Maxwellian atomic data are automatically propagated into kappa distribution 
data obtained using the decomposition approach. This is in contrast to reverse-engineering 
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methods, for which it is necessary to recalculate the relevant rate coefficients. Additionally, 
the uncertainties with the decomposition approach mainly reflect the accuracy with which 
the Maxwellian sum approximates the kappa distribution. Thus, hnding a more accurate 
approximation will improve all of the data. Finally, any value of k can, in principal, be 
approximated to the same accuracy so that the atomic data derived for small k can be as 
accurate as that derived for large k. 

We have applied the Maxwellian decomposition method to several examples, mainly 
focusing on the CSD. First, we incorporated EIMI into CIE calcualations and found that 
EIMI can signihcantly change the equilibrium CSD. The effect is greatest for small values 
of K, which makes sense as such distributions have the largest nonthermal tails. Next, we 
studied the evolution of the CSD in a plasma in which both k and can vary. For small 
values of k, the CSD is predicted to evolve towards equilibrium faster than those for larger 
values of k or Maxwellian plasmas. The time dependent CSD equation was solved for 
several cases that model flares as generating a kappa EED and the time evolution of the 
various charge states was found. The charge states in a kappa distribution plasma evolve 
towards a different equilibrium than those for a thermal distribution. EIMI was found to 
speed up the rate at which the plasma ionizes up through the lower charge states. 

The method of approximating kappa distributions as a sum of Maxwellians is not new, 
but seems to have fallen out of favor recently in comparison to other methods. This neglect 
is not justihed. The accuracy, ease of use, and portability of the method make it a useful 
tool for spectroscopic analysis. The same methods can also be applied to other nonthermal 
distributions beyond kappa distributions. The main issue is to hud systematic methods for 
performing the Maxwellian decomposition to high accuracy. 

This work was supported in part by the NASA Living with a Star program grant 


NNX15AB71G. 
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Table 1. Fitting Parameters for Maxwellian Decomposition of Kappa Distributions with 

K = 1.7, 2, and 3 (Equation [13]) 


j 

K, = 

aj 

1.7 

Cj 

aj 

= 2 

D 

K 

Qj 

= 3 

Cj 

0 

5.00(-2) 

9.5183(-2) 

1.02(-1) 

5.9896(-2) 

1.16(-1) 

2.5922(-3) 

1 

1.24(-1) 

3.5911(-1) 

1.68(-1) 

1.3330(-1) 

2.56(-l) 

9.3798(-2) 

2 

3.27(-l) 

3.2983(-l) 

3.08(-l) 

3.3109(-1) 

3.54(-l) 

1.0364(-1) 

3 

9.23(-l) 

1.4888(-1) 

4.41(-1) 

5.2666(-2) 

5.51(-1) 

3.3185(-1) 

4 

2.798(0) 

5.4011(-2) 

7.73(-l) 

2.7887(-l) 

8.72(-l) 

1.8873(-1) 

5 

4.266(0) 

-9.8173(-3) 

2.043(0) 

1.0852(-1) 

1.391(0) 

1.9666(-1) 

6 

7.695(0) 

2.0397(-2) 

5.568(0) 

2.5088(-2) 

2.920(0) 

6.8733(-2) 

7 

1.320(1) 

-8.5302(-3) 

8.611(0) 

4.1911(-3) 

6.437(0) 

1.1027(-2) 

8 

2.0017(1) 

1.6425(-2) 

2.2182(1) 

1.3765(-2) 

1.1726(1) 

2.4131(-3) 

9 

2.4556(1) 

-8.5330(-3) 

2.6320(1) 

-1.1296(-2) 

2.4583(1) 

3.7881(-4) 

10 

4.5069(1) 

2.4445(-3) 

3.8883(1) 

3.5490(-3) 

4.1266(1) 

1.9485(-4) 

11 

2.22953(2) 

1.1895(-3) 

2.08310(2) 

5.0451(-4) 

1.88810(2) 

-4.1575(-5) 

12 

3.52864(2) 

-1.4590(-3) 

3.11017(2) 

-2.2212(-4) 

7.95476(2) 

2.9298(-5) 

13 

5.32171(2) 

7.6344(-4) 

9.30643(2) 

-1.0136(-4) 



14 

1.27669(3) 

3.0240(-4) 

2.45618(3) 

1.7277(-4) 



15 

3.84582(3) 

-1.9796(-4) 






Note. — The entries of the form a(b) are shorthand for a x 10^. 
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Table 2. Same as Tabled] but for k = 4, 5, and 7 


K = 4 K = 5 K = 7 


J 

ttj 


ttj 

Cj 

Qj 


0 

3.29(-l) 

1.1589(-1) 

4.62(-l) 

2.6404(-l) 

4.92(-l) 

1.6817(-1) 

1 

5.98(-l) 

3.9529(-l) 

8.64(-l) 

4.9170(-1) 

8.12(-1) 

4.8515(-1) 

2 

1.053(0) 

3.2316(-1) 

1.577(0) 

2.0295(-l) 

1.320(0) 

2.8210(-1) 

3 

1.844(0) 

1.2715(-1) 

2.921(0) 

3.7629(-2) 

2.195(0) 

5.9684(-2) 

4 

3.302(0) 

3.2476(-2) 

5.936(0) 

3.5942(-3) 

3.989(0) 

4.8443(-3) 

5 

6.411(0) 

5.5656(-3) 

2.0618(1) 

8.1131(-5) 

1.1406(1) 

5.3502(-5) 

6 

1.5520(1) 

4.7317(-4) 






Table 3. Same as Tabled] but for k, = 10, 15, and 20 


j 

K 

Qj 

= 10 

Cj 

K 

ttj 

= 15 

K ■ 

ttj 

= 20 

Cj 

0 

5.18(-1) 

9.5556(-2) 

6.37(-l) 

1.9206(-1) 

6.06(-l) 

5.9165(-2) 

1 

7.84(-l) 

4.5424(-l) 

9.61(-1) 

6.3165(-1) 

8.30(-l) 

4.4999(-l) 

2 

1.190(0) 

3.7406(-l) 

1.526(0) 

1.7558(-1) 

1.139(0) 

4.2849(-l) 

3 

1.920(0) 

7.4674(-2) 

3.966(0) 

7.0226(-4) 

1.643(0) 

6.2175(-2) 

4 

5.058(0) 

4.8765(-3) 

5.633(0) 

1.3855(-5) 

3.475(0) 

1.7974(-4) 

5 

5.550(0) 

-3.3968(-3) 

7.0359(1) 

-8.8205(-6) 
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Table 4. Same as Table [H but for k = 25, 30, and 33 


j 

K 

ttj 

= 25 

K 

aj 

II 

CO 

o 

K 

aj 

= 33 

0 

6.62(-l) 

1.0281(-1) 

7.29(-l) 

1.8056(-1) 

7.16(-1) 

1.2455(-1) 

1 

9.26(-l) 

6.486(-l) 

9.76(-l) 

6.3110(-1) 

9.39(-l) 

5.8530(-l) 

2 

1.328(0) 

2.4775(-l) 

1.333(0) 

1.8601(-1) 

1.185(0) 

2.1208(-1) 

3 

2.827(0) 

8.4738(-4) 

1.873(0) 

2.2955(-3) 

1.406(0) 

7.8024(-2) 

4 

1.4319(1) 

-1.5625(-5) 

8.438(0) 

5.3631(-5) 

6.374(0) 

5.0427(-5) 

5 

6.9158(1) 

3.4698(-6) 

3.145(1) 

-1.8321(-5) 

8.7117(1) 

-7.4012(-6) 


Table 5. Same as Table [U but for k = 50 and 100 


j 

tl = 

aj 

50 

K = 

Qj 

100 

Cj 

0 

7.23(-l) 

5.0910(-2) 

8.43(-l) 

2.1409(-1) 

1 

9.16(-1) 

6.0722(-l) 

1.023(0) 

7.2273(-l) 

2 

1.190(0) 

3.4184(-1) 

1.273(0) 

6.3198(-2) 

3 

1.2300(1) 

4.6799(-5) 

2.0710(1) 

-2.1750(-5) 

4 

1.00200(2) 

-1.6831(-5) 
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Table 6. Accuracies of Fitting Parameters in Tables HHS] 


K 

Amax/^BT'K’^ 

Amax for E < Ajiiax 

Ej\cj\ 

1.7 

632 

0.0214 

1.057 

2 

330 

0.0250 

1.023 

3 

77.7 

0.0194 

1.000 

4 

41.5 

0.0020 

1.000 

5 

29.7 

0.0032 

1.000 

7 

21.1 

0.0009 

1.000 

10 

16.7 

0.0032 

1.007 

15 

14.2 

0.0262 

1.000 

20 

13.1 

0.0017 

1.000 

25 

12.5 

0.0165 

1.000 

30 

12.2 

0.0077 

1.000 

33 

12.0 

0.0055 

1.000 

50 

11.5 

0.0282 

1.000 

100 

11.0 

0.0020 

1.000 


is the energy below which 99.99% of the 
particles are found, i.e., /^(A; k, T^) dA = 


0.9999. 






Table 7. Fitting Kappa Distributions for Continuously Varying k = 1.7 - 100 


Range 

aj 

do 

di 

d2 

do 

6^4 ^5 (Iq 

1.7 < K < 2.0 

0.049 

1.56729(0) 

-1.59895(0) 

4.07812(-1) 



n 

0.080 

8.35027(0) 

-1.08223(1) 

4.68258(0) 

-6.76624(-l) 


n 

0.168 

-3.94995(-l) 

2.00317(0) 

-1.33084(0) 

2.42976(-l) 


n 

0.384 

-7.98933(0) 

1.06867(1) 

-4.51282(0) 

6.33092(-l) 


n 

0.944 

1.03347(1) 

-1.42760(1) 

6.50682(0) 

-9.48923(-l) 


n 

1.118 

-3.37385(1) 

5.75809(1) 

-3.66138(1) 

1.03267(1) 

-1.09760(0) 

II 

2.267 

2.54497(0) 

-3.50125(0) 

1.59975(0) 

-2.34561(-1) 


II 

4.362 

-1.92314(0) 

2.58894(0) 

-1.12945(0) 

1.61818(-1) 


II 

12.253 

3.98211(-1) 

-5.25104(-1) 

2.35482(-l) 

-3.49376(-2) 


II 

38.992 

-3.94167(0) 

7.29473(0) 

-5.04540(0) 

1.54697(0) 

-1.77514(-1) 

II 

153.710 

1.56669(0) 

-1.54177(0) 

3.79314(-1) 



II 

199.674 

-2.03669(0) 

1.99049(0) 

-4.86039(-l) 



II 

393.810 

7.76167(-1) 

-7.39125(-1) 

1.75439(-1) 



II 

997.366 

-1.76300(-1) 

1.61646(-1) 

-3.67022(-2) 



2.0 < K < 2.3 

0.049 

-4.77764(-2) 

4.21978(-2) 

-9.00028(-3) 



II 

0.080 

8.35027(0) 

-1.08223(1) 

4.68258(0) 

-6.76624(-l) 


II 

0.168 

-3.94995(-l) 

2.00317(0) 

-1.33084(0) 

2.42976(-l) 


II 

0.384 

-7.98933(0) 

1.06867(1) 

-4.51282(0) 

6.33092(-l) 


II 

0.944 

1.03347(1) 

-1.42760(1) 

6.50682(0) 

-9.48923(-l) 


II 

1.118 

-3.37385(1) 

5.75809(1) 

-3.66138(1) 

1.03267(1) 

-1.09760(0) 

II 

2.267 

2.54497(0) 

-3.50125(0) 

1.59975(0) 

-2.34561(-1) 


II 

4.362 

-1.92314(0) 

2.58894(0) 

-1.12945(0) 

1.61818(-1) 


II 

12.253 

3.98211(-1) 

-5.25104(-1) 

2.35482(-l) 

-3.49376(-2) 


II 

38.992 

-3.94167(0) 

7.29473(0) 

-5.04540(0) 

1.54697(0) 

-0.177514(-1) 

II 

153.710 

3.13273(-2) 

-2.96564(-2) 

7.09938(-3) 



2.3 < K < 2.4 

0.049 

-4.77764(-2) 

4.21978(-2) 

-9.00028(-3) 



II 

0.080 

8.35027(0) 

-1.08223(1) 

4.68258(0) 

-6.76624(-l) 


II 

0.168 

-3.94995(-l) 

2.00317(0) 

-1.33084(0) 

2.42976(-l) 


II 

0.384 

-7.98933(0) 

1.06867(1) 

-4.51282(0) 

6.33092(-l) 


II 

0.944 

1.03347(1) 

-1.42760(1) 

6.50682(0) 

-9.48923(-l) 


II 

1.118 

-3.37385(1) 

5.75809(1) 

-3.66138(1) 

1.03267(1) 

-1.09760(0) 








Table 7—Continued 


Range 

aj 

do 

di 

d2 

ds 

^4 

d?5 


„ 

2.267 

2.54497(0) 

-3.50125(0) 

1.59975(0) 

-2.34561(-1) 




" 

4.362 

-1.92314(0) 

2.58894(0) 

-1.12945(0) 

1.61818(-1) 




" 

12.253 

3.98211(-1) 

-5.25104(-1) 

2.35482(-l) 

-3.49376(-2) 




" 

38.992 

-3.94167(0) 

7.29473(0) 

-5.04540(0) 

1.54697(0) 

-1.77514(-1) 


2.4 < 

K < 3.7 

0.116 

3.96920(0) 

-5.61929(0) 

3.19670(0) 

-9.11843(-1) 

1.30249(-1) 

-7.44621(-3) 


" 

0.256 

5.64248(0) 

-5.33407(0) 

1.93157(0) 

-3.12193(-1) 

1.78935(-2) 

2.02901(-4) 


" 

0.354 

-8.32819(0) 

9.51261(0) 

-4.18395(0) 

8.93371(-1) 

-9.16849(-2) 

3.51840(-3) 


" 

0.551 

3.37016(0) 

-2.76432(0) 

7.98392(-l) 

-3.48088(-2) 

-1.90868(-2) 

2.28428(-3) 


" 

0.872 

-3.64079(0) 

3.60739(0) 

-1.27424(0) 

1.93318(-1) 

-7.60047(-3) 

-5.29236(-4) 


" 

1.391 

9.04059(-l) 

-6.29854(-l) 

1.24991(-1) 

2.96435(-2) 

-1.25896(-2) 

1.13885(-3) 


" 

2.920 

-1.04513(0) 

1.22861(0) 

-5.33511(-1) 

1.14042 (-1) 

-1.19138(-2) 

4.74984(-4) 


" 

6.437 

3.96190(-1) 

-3.53759(-l) 

1.23179(-1) 

-1.88395(-2) 

7.98385(-4) 

4.81100(-5) 


" 

11.726 

-4.96928(-l) 

6.30421(-1) 

-3.20607(-l) 

8.18266(-2) 

-1.04245(-2) 

5.28095(-4) 


" 

24.583 

2.28994(-l) 

-2.77800(-l) 

1.37518(-1) 

-3.45280(-2) 

4.36068(-3) 

-2.20187(-4) 

3.7 < 

K < 4.5 

0.302 

1.96013(0) 

-1.10591(0) 

2.18455(-1) 

-1.48857(-2) 




" 

0.561 

-2.15731(0) 

1.69698(0) 

-3.70337(-l) 

2.65640(-2) 




" 

1.043 

3.80733(0) 

-2.45503(0) 

5.63680(-l) 

-4.15712(-2) 




" 

1.590 

-4.73982(0) 

3.19782(0) 

-7.01062(-1) 

5.08460(-2) 




" 

2.582 

2.76133(0) 

-1.77615(0) 

3.89392(-l) 

-2.83457(-2) 




" 

6.215 

-1.73113(0) 

1.17320(0) 

-2.60640(-l) 

1.90386(-2) 




" 

9.626 

1.43348(0) 

-9.54201(-1) 

2.09719(-1) 

-1.52106(-2) 




" 

17.292 

3.34049(-l) 

2.23321(-1) 

-4.91735(-2) 

3.56517(-3) 



4.5 < 

K < 5.2 

0.302 

6.04673(0) 

-5.53465 (0) 

2.11477(0) 

-4.14777(-1) 

4.136333(-2) 

-1.66875(-3) 


" 

0.561 

-5.59121(0) 

4.95522(0) 

-1.52890(0) 

2.09511(-1) 

-1.08252(-2) 



" 

1.043 

1.34115(1) 

-1.13733(1) 

3.66504(0) 

-5.20241(-1) 

2.76641(-2) 



" 

1.590 

-2.08424(1) 

1.79968(1) 

-5.79182(0) 

8.2769(-l) 

-4.43552(-2) 



" 

2.528 

1.72935(1) 

-1.49928(1) 

4.88582(0) 

-7.06458(-l) 

3.82436(-2) 



" 

6.215 

-2.53090(1) 

2.23218(1) 

-7.35343(0) 

1.07306(0) 

-5.85483(-2) 


5.2 < 

K < 7.3 

0.397 

1.69355(0) 

-7.38559(-l) 

1.29662(-1) 

-1.05755(-2) 

3.32008(-4) 



" 

0.662 

-1.64244(0) 

1.04587(0) 

-1.98344(-1) 

1.67171(-2) 

-5.32776(-4) 



" 

1.130 

1.36114(0) 

-6.07907(-l) 

1.29423 (-1) 

-1.12748(-2) 

3.60750(-4) 



.. 

2.133 

-7.80854(-l) 

4.69366(-l) 

-8.92911(-2) 

7.22596(-3) 

-2.15370(-4) 







Table 7—Continued 


Range 

aj 

do 

di 

C?2 

do 

6^4 

C?5 

do 


n 

6.708 

-8.22231(-1) 

8.95811(-1) 

-3.17094(-1) 

5.10962(-2) 

-3.93547(-3) 

1.19768(-4) 



n 

10.978 

7.02190(-1) 

-6.65658(-l) 

1.87408(-1) 

-2.06380(-2) 

7.72861(-4) 




n 

19.325 

-3.57098(0) 

2.81704(0) 

-8.50939(-l) 

1.25010(-1) 

-9.07755(-2) 

2.65483(-4) 



■■ 

49.952 

4.506(-l) 

-3.080(-l) 

7.376(-2) 

-7.380(-3) 

2.596(-4) 



7.3 < 

A 

00 

0.397 

1.69355(0) 

-7.38559(-l) 

1.29662(-1) 

-1.05755(-2) 

3.32008(-4) 




■■ 

0.662 

-1.64244(0) 

1.04587(0) 

-1.98344(-1) 

1.67171(-2) 

-5.32776(-4) 




■■ 

1.130 

1.36114(0) 

-6.07907(-l) 

1.29423(-1) 

-1.12748(-2) 

3.60750(-4) 




■■ 

2.133 

-7.80854(-l) 

4.69366(-l) 

-8.92911(-2) 

7.22596(-3) 

-2.15370(-4) 




■■ 

6.708 

1.18523(-1) 

-4.04953(-2) 

4.58953(-3) 

-1.68968(-4) 





■■ 

10.978 

-4.96980(-2) 

1.78065(-2) 

-2.08221(-3) 

7.77779(-5) 




VI 

00 

K < 11.5 

0.452 

1.30350(0) 

-4.49713(-1) 

6.88978(-2) 

-5.82877(-3) 

2.81686(-4) 

-7.30064(-6) 

7.88427(-8) 


■■ 

0.694 

-1.66061(0) 

8.73637(-l) 

-1.52021(-1) 

1.38576(-2) 

-7.05050(-4) 

1.89861(-5) 

-2.11219(-7) 


■■ 

1.074 

2.24466(0) 

-9.18802(-1) 

1.76499(-1) 

-1.69900(-2) 

8.96474(-4) 

-2.47898(-5) 

2.81438(-7) 


■■ 

1.791 

-1.56419(0) 

7.86300(-l) 

-1.45517(-1) 

1.38959(-2) 

-7.32919(-4) 

2.03092(-5) 

-2.31248(-7) 


■■ 

5.315 

4.89878(-l) 

-1.58942(-1) 

2.12165(-2) 

-1.534327(-3) 

6.24996(-5) 

-1.10994(-6) 



■■ 

9.741 

-3.01963(-1) 

6.00271(-2) 

-2.08105(-3) 

-1.66306(-4) 

7.59614(-6) 




■■ 

22.641 

1.01885(-1) 

1.06359(-2) 

-7.32292(-3) 

7.19906(-4) 

-1.96401(-5) 




■■ 

74.007 

-4.89289(-3) 

-1.91647(-2) 

5.18355(-3) 

-4.33200(-4) 

1.11707(-5) 



11.5 < 

K < 13.5 

0.452 

1.30350(0) 

-4.49713(-1) 

6.88978(-2) 

-5.82877(-3) 

2.81686(-4) 

-7.30064(-6) 

7.88427(-8) 


■■ 

0.694 

-1.66061(0) 

8.73637(-l) 

-1.52021(-1) 

1.38576(-2) 

-7.05050(-4) 

1.89861(-5) 

-2.11219(-7) 


■■ 

1.074 

2.24466(0) 

-9.18802(-1) 

1.76499(-1) 

-1.69900(-2) 

8.96474(-4) 

-2.47898(-5) 

2.81438(-7) 


■■ 

1.791 

-1.56419(0) 

7.86300(-l) 

-1.45517(-1) 

1.38959(-2) 

-7.32919(-4) 

2.03092(-5) 

-2.31248(-7) 


■■ 

5.315 

2.50719(-2) 

-4.96470(-3) 

3.31831(-4) 

-7.42211(-6) 




13.5 < 

A 

00 

b 

0.452 

1.30350(0) 

-4.49713(-1) 

6.88978(-2) 

-5.82877(-3) 

2.81686(-4) 

-7.30064(-6) 

7.88427(-8) 


■■ 

0.694 

-1.66061(0) 

8.73637(-l) 

-1.52021(-1) 

1.38576(-2) 

-7.05050(-4) 

1.89861(-5) 

-2.11219(-7) 


■■ 

1.074 

2.24466(0) 

-9.18802(-1) 

1.76499(-1) 

-1.69900(-2) 

8.96474(-4) 

-2.47898(-5) 

2.81438(-7) 


■■ 

1.791 

-1.56419(0) 

7.86300(-l) 

-1.45517(-1) 

1.38959(-2) 

-7.32919(-4) 

2.03092(-5) 

-2.31248(-7) 

18.0 < K 

0.411 

4.7555(-2) 

-4.1569(-3) 

1.5348(-4) 

-3.0868(-6) 

3.4844(-8) 

-2.0732(-10) 

5.0592(-13) 


■■ 

0.755 

7.6317(-1) 

-3.9928(-2) 

1.2674(-3) 

-2.3391(-5) 

2.4994(-7) 

-1.4319(-9) 

3.998(-12) 


■■ 

1.029 

-2.0600(-l) 

6.9082(-2) 

-2.2280(-3) 

4.1540(-5) 

-4.4691(-7) 

2.5731(-9) 

-6.1319(-12) 


" 

1.551 

3.9528(-l) 

-2.4997(-2) 

8.0719(-4) 

-1.5062(-5) 

1.62134(-7) 

-9.3385(-10) 

2.2261(-12) 
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Note. — Here Cj{K) = dni^^- The entries of the form a(b) are shorthand for a x 10^. 
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Fig. 1.— The maximum absolute relative error i?max versus the sum of the magnitudes of 
the Cj, kil- A random number generator was used to choose the set of aj for n = 7. This 
shows that there is a strong positive correlation between these quantities, indicating that 
minimizing \cj\ tends to also minimize the error. 
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Fig. 2.— Example of a Maxwellian decomposition fit to /k(E) for k = 2 and = 10 (the 
units are arbitrary). The dashed curve is the kappa distribution, while the solid curve shows 
the result obtained by summing Maxwellians. The vertical dotted line indicates the energy 
at which 99.99% of the particles have a smaller energy. This corresponds to 330 ksTi^ for 
n = 2. The lower panel shows i?err, which is at most 0.025 for E < 330 k-QT^. 
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Fig. 3.— The maximum relative error i^max between the kappa distribution and the sum 
of Maxwellians approximated using the values in Table [71 The solid line indicates the i?max 
found in the energy range from zero up to the energy that contains 99.9% of the kappa 
distribution and the dotted line shows i?max for the higher energy that contains 99.99% of 
the particles. 
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104 jq5 jq6 jq7 jq8 


T. (K) 


Fig. 4.— The ionization rate coefficient ai for Fe^"*" forming Fe^"*" as a function of for 


K = 2. The solid curve shows our result using the Maxwellian deco mposition method, whi 


the dashed line indicates the result using the KAPPA package of 


Dzifcakova et ah 


(120151). 


The lower panel indicates the relative error given by (KAPPA result - our result) / our 
result. 
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Fig. 5.— Same as Figure 01 but for the total recombination rate coefficient or for Fe^"*" 
forming Fe^^. The rate coefficient includes both RR and DR. 
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Fig. 6.— The CSD for Fe with k = 2 as a function of log[T«;(K)]. The top panel shows 
the relative abundances of the various charge states, with our results indicated by the solid 


curves and the abundances given by the KAPPA packag e shown as dashec 


panel shows the ratio of our (New) results to those of 


Dzifcakova et af 


curves . The lower 
fania Old). The 


curves are plotted only for where the abundances are greater than 1%. 
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Fig. 7.— The intensity ratio of two pairs of O VI lines as a function of k. The emissivities 
were calculated for = 1 MK and Ug = 10® cni“^. The solid curves show our results, 
which were calculated in a straightforward way using built-in CHIANTI functions (see text 
for details). The filled circles indicate the ratios obtained using the KAPPA package. 
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Fig. 8.— The CSD for Fe as a function of log [Tk(K)] for k, = 25. The top panel shows 
the relative abundances of the varous charge states. The solid curve indicates the results 
when EIMI is included in the calculation and the dashed curve illustrates the CSD when 
only single ionization is considered. The lower panel shows the ratio of the abundances 
including (EIMI) or neglecting EIMI (Single). The curves are plotted only for where the 
abundances are greater than 1%. 
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Fig. 9.— Same as Figured but for n 


1.7. 
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Fig. 10.— The equilibration timescale of the CSD as a function of k at several values of T^. 
This shows that the CSD evolves more rapidly for smaller values of k, which are the most 
non-thermal. 
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Fi g. 11.— The evolution of the CSD in response to a flare, based on parameters inferred 


by 


Kasparova fc Karlickvi (1200911 from hard X-ray observations. The initial condition is a 


Maxwellian plasma at 1 MK. The solid curve shows the calculated abundances, including 
EIMI, for a flare that produces electrons with a distribution described by k = 8 and = 
15 MK. The dashed curve illustrates the same conditions, but neglecting EIMI processes. 
The dotted curve illustrates the results if the plasma were to remain Maxwellian and only 
have an increased temperature of 15 MK. 
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